PARTIAL DERIVATIVE APPROACH FOR OPTION PRICING 
IN A SIMPLE STOCHASTIC VOLATILITY MODEL 



MIQUEL MONTERO 
Departament de Fisica Fonamental, Universitat de Barcelona, 
Diagonal 647, 08028 Barcelona, Spain 
E-mail: miquel . montero @ ub . edu 



Abstract 

We study a market model in which the volatility of the stock may 
jump at a random time r from a fixed value a a to another fixed value 
<7j,. This model was already described in the literature. We present a new 
approach to the problem, based on partial derivative equations, which 
gives a different perspective to the problem. Within our framework we 
can easily consider several prescriptions for the market price of volatility 
risk, and interpret their financial meaning. Thus, we recover solutions 
previously cited in the literature as well as obtain new ones. 



1 Introduction 

The problem of pricing financial derivatives was already present in the aim 
of the early works in Mathematical Finance. Bachelier in 1900 proposed the 
arithmetic Brownian motion for the dynamical evolution of stock prices as a 
first step towards obtaining a price for options [3]. Nevertheless the interest 
on this problem has increased remarkably in the past twenty years, after the 
publication of the works of Black and Scholes [2], and Merton [12]. The Black- 
Scholcs model has been broadly used by practitioners thereafter, mainly due 
to its mathematical simplicity. It is well established, however, that this model 
fails to explain some statistical features shown in real markets. In particular, 
there are solid evidences pointing to the necessity of relaxing the assumption, 
present in the Black-Scholes model, that a constant volatility parameter drives 
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the stock price. One of the tests more commonly used is based on a conceptually 
simple principle. Since the Black-Scholes price is a monotonous function on 
its arguments, the formula can be inverted in order to compute the implied 
volatility, the volatility that will reproduce the actual market conditions. The 
usual result is that the implied volatility is not constant, but a U-shaped function 
of the moneyness, whose minimum is at moneyness near to one — i.e. when the 
current price of the underlying is close to the strike. This departure from the 
Black-Scholes model is known as the smile effect, and it is well documented in 
the literature [5]. 

Many models have been developed with the purpose of avoiding this inac- 
curate feature. We will mention here only a few of them. Merton itself [13] 
proposed a model in which the volatility was a deterministic function of time. 
Cox and Ross [4] presented some alternative proposals that can be thought 
as models in which the volatility is stock-dependant. These and other similar 
contributions lead to a framework in which all the option risk comes from the 
fluctuations in the price of the underlying. In practical situations, however, it 
seems that this description is not sophisticated enough for explaining the actual 
changes in the level of volatility. Some authors have then suggested that the 
evolution of the volatility is driven by its own stochastic equation. Among these 
models of stochastic volatility we find works that are historically noteworthy: 
Hull and White [9] proposed a model where the squared volatility also follows 
a log-normal diffusion equation, independent of the stock price. Wiggins [17] 
extended this idea and considered that the underlying and the volatility con- 
stitute a two-dimensional system of correlated log-normal random processes. 
Scott [15], but specially Stein and Stein [16] assumed that the instantaneous 
volatility follows a random mean-reverting process: an independent arithmetic 
Ornstein-Uhlcnbcck process. Masoliver and Perello [10] relaxed this assumption, 
and introduced correlation in the two-dimensional Wiener process. Heston [8] 
turned the arithmetic model into a square-root correlated process. 

All these seminal papers have in common that they model the stochastic be- 
haviour of the volatility as a diffusion process. Naik [14] developed a model in 
which the volatility can have only two known values, and the market switches 
back and forth between them, in a random way. This set-up can be used to 
model a market with high and low volatility periods. Herzel [7] studied a sim- 
plified version of this problem, in which the volatility, at the most, can jump 
once. This is a suitable model for encoding a market that may undergo a severe 
change in volatility only if some forthcoming event takes place. Since options 
have a limited lifetime, this seems not to be a very restrictive limitation. Herzel 
solved the problem of pricing the options using probability arguments, and 
showed that his model can account for the smile effect. 

We present here a different approach for obtaining fair option prices under 
Hcrzcl's conditions. We will employ a technique of broad use both in research 
papers (e.g. in [8]) and reference books (e.g. in [18]) on this topic: we determine 
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the partial derivative equations that the option price must fulfil, according to the 
ltd convention, and solve it with the appropriate constrains. This scheme eases 
the task of considering several prescriptions for the market price of volatility 
risk, and leads to a plain way of interpreting the financial meaning of each of 
them. 

The paper is structured as follows: in Sec. 2 we present the general market 
model and specify the differential equations that govern the traded securities. 
In Sec. 3 we study the way of obtaining a complete market. In Sec. 4 we explore 
the consequences of demanding that the market admits no arbitrage. In Sec. 5 
we present explicit solutions for different market prices of the volatility risk. 
Section 6 contains actual numerical examples of these solutions and a financial 
interpretation of the results. The conclusions are drawn in Sec. 7. The paper 
ends with Appendix A, where we detail the way we have followed for finding of 
one of the new solutions we have introduced. 

2 The Market Model 

Let us begin with the general description of our set-up. We will assume that in 
our market there is at least a non-deterministic traded stock, S. The evolution of 
the price of this stock, from So at t — t , is governed by the following differential 
equation" : 

dS , . _ 

— = ndt + adW, 
o 

where W(t — to) is a one dimensional Brownian motion, with zero mean and 
variance equal to t — t 0} \x is a constant parameter, and a, the volatility, is a 
stochastic quantity. The model assumes that the volatility have initially a given 
value (T , and that at most it may change to a different value ct at instant 
t > to: 

a(t;r) = a a \ t<T + cr b l t > T = a a + (a b - a a )l t > T , (1) 

where denotes the indicator function, which assigns the value 1 to a true 
statement, and the value to a false statement. The time r in which such 
transition occurs is random and we will assume that it follows an exponential 
law: 

P(to<T<t) = l-e- x{t - to \ 

Note that with the previous definition, A is just the inverse of the mean transition 
time, E[t] = A -1 . 

We are also assuming that we will be capable of concluding whether the 
transition has taken place or not. This assumption does not imply that we can 
directly measure the value of cr, but that there exists a way to determine if 
t > t. This can be easily understood from the point of view of a practitioner. 

"Throughout our exposition we will not specify the explicit dependence of the involved 
magnitudes, except if this may lead to confusion. 
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Let us suppose, for instance, that we are expecting that a relevant financial 
announcement is done. We do not know for sure when it will happen, but we 
bclivc that this new information will affect the level of volatility in our market. 
Even though we may not perform an instantaneous measure of the volatility in 
order to check the actual effect of the news, if they are published we will know 
it. We will return to this issue later. 

Upon the underlying stock S, we will define a new traded asset: the option 
C. The price of this option will depend explicitly on the moment to in which 
we will decide to evaluate it, on the current stock price So and on the level of 
volatility <7o, but also on a set of peculiar parameters which we will label with 
a single symbol, n. These parameters are the contract specifications that will 
characterize the option: the maturity time or the striking price, among other 
possibilities. This framework covers the European put and call options, e.g. 
the vanilla options or the binary options; also the American options, but does 
not include more exotic derivatives, such as the Asian options or the lookback 
options. 

The differential of the option price C = C(t, S, a; k) has, according to the 
Ito convention, the following expression: 

1 AC 

dC = dtCdt + d s CdS + -a 2 S 2 d 2 ss Cdt + da, (2) 

2 crh — a a 

where 

AC = C(t, S, a b ; k) - C(t, S, a a ; k). 

The last term in Eq. (2) condenses the innovation with respect to the classical 
Black-Scholes expression, and represents the contribution of the randomness in 
the volatility to the dynamics of the option price. Note that this extra term is 
a product of the finite difference version of the derivative of C with respect to 
a, and da. In order to obtain an alternative expression for this object, we will 
simply differentiate Eq. (1): 

da = (a b - a a )dl t > T . (3) 

The differential of a indicator may seem a bizarre object. However, it is math- 
ematically well defined, as we will shortly show. We can decompose this differ- 
ential in two terms: 

dl t > T = Xl t<T dt - XdG. (4) 
The first term is regular, and the second involves function G, 

G = tl t<T + (j-^j !t>T, 

which is proven to be a right continuous with left limits martingale. In order 
to ease the notation, nevertheless, we will keep the differential of the indicator 



PDE for a Stochastic Volatility Model 



5 



in its early form, and use the referred decomposition only when it can clarify 
the problem. We must stress however that rfl t > T is a stochastic magnitude, 
independent of dW . Since da does not directly contribute on the variation of 
the stock price dS, we can foreseen that there is a source of risk that cannot 
be explained in terms of the random evolution of the underlying asset. We 
postpone nonetheless the discussion of this issue, since it will be the matter of 
the next Section. 

Before that we want to point out that there is also a third kind of security 
traded in the market, a free-risk monetary asset B, which satisfies the corre- 
sponding equation: 

dB = rBdt. (5) 

This security will allow us to borrow money when we need it, and it will provide 
a secure resort in order to keep the excess of cash if that is the case. In particular, 
it makes possible both the self-financing strategy, which allows closed portfolios, 
and the net-zero investment, the composition of a portfolio with no net value. 

3 Completeness of the Market 

Let us face the problem of the completeness of the market. It is notorious that 
the market will be complete if we can construct for every security the so-called 
replicating portfolio, i.e. a portfolio that mimics the behaviour of the asset. 
We have argued in the previous Section that not all the influence of a in the 
price of the option can be explained through S. We need then another security 
that can account for this component of the global risk. Instead of introducing a 
new traded asset depending only on da, with no clear financial interpretation, 
we have decided to use a secondary option D(t, S, a; n'): a derivative of the 
same nature of C(t, S,a; n), but with a different set of contract specifications. 
This add-on completes the market if we are allowed to borrow money at a fixed 
interest rate whenever we need it, or to buy zero-coupon bonds in the case we 
obtain a surplus of cash. Thus we can write down C as a combination of 5 
shares S, </> units of the riskless security B, and tp secondary options D: 

C = 5S + (j)B + ipD. 

The variation in the value of both portfolios fulfills 

dC = SdS + 4>dB + ipdD, 

where we have taken into account two capital facts. On one hand S, <p and ip 
are nonanticipating functions of S and D, e.g. dS, d<p and dip do not depend 
on the new random information in dW and da. On the other hand, we adopt 
a self-financing strategy, in which there is no net cash flow entering or leaving 
the replicating portfolio [6]: 

SdS + Bd<t> + DdiP = 0. 
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We will replace dC by the expression in (2), and we will take into account 
the properties shown in (3) and (5), to finally obtain: 

dtCdt + d s CdS + ^cr 2 S 2 d 2 ss Cdt + ACdl t > T = SdS + r^Bdt + ifjdD. (6) 
We can proceed with dD in an analogous way 

dD = d t Ddt + d s DdS + ^a 2 S 2 d 2 ss Ddt + ADdl t > T , (7) 
where the natural definition of AD, 

AD = D{t, S, <r h ; k') - D(t, S, a a ; n'), 

has been used. In order to recover a deterministic partial differential equation 
we must guarantee that all the terms containing the stochastic magnitudes dS 
and dl t > T mutually cancel out. Thus we must demand that 

8 S C = 6 + yjd s D, 

condition named delta hedging, and also that 

AC = ijAD, 

which is usually referred as vega hedging, or sometimes as psi hedging [11]. 
The previous hedging conditions reduce Eq. (6) to 

d t Cdt + \a 2 S 2 d 2 ss Cdt = rcpBdt + ^f- (d t Ddt + \o 2 S 2 d 2 ss Ddt) , (8) 

expression that still involves B, which is not an inner variable of the option 
prices C and D. This problem can be fixed using together the definition of the 
portfolio and the psi hedging, 

^B = C-5S-^D = C- (d s C - ^d s D^j S - |^D. 

Thus, the replacement of cf)B in Eq. (8) leads to 

d t C+^a 2 S 2 d 2 ss C-rC+rSd s C = ^ (d t D + ^a 2 S 2 d 2 ss D - rD + rSd s D^j . 

This formula implies the existence of an arbitrary function \ = x(A S, cr), which 
uncouples the problem of finding C and D: 

X (t, S,a) = ^ (d t C + \cr 2 S 2 d 2 s C -rC + rSdsC^j . (9) 

Obviously the same formula, just replacing C by D, is valid for the secondary 
option. This fact proves that the option D completes the market indeed [1, 11]. 
Note that if \{t, S,<r) — we recover the classical Black-Scholes equation. 
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4 Arbitrage-free scenario 

We have to state some criterion before we can choose a valid candidate for 
X(t, S, a). We will first determine the meaning of this arbitrary function. Con- 
sider Y, a portfolio which involves shares, bonds, one primary option, and sec- 
ondary options: 

Y = C + 8S + 4>B + $D. 

Moreover, the relative amount of each security is such that, at the beginning, 
the portfolio has no net value, i.e. Y = 0. In fact, we will also demand that 
the changes in the value of the portfolio arc not coming from a cash flow. The 
nonanticipating nature of 8, <j>, and makes that dY takes the following form: 

dY = dC + 8dS -r(C + SS + ^L>) dt. 

Finally, we can use Eqs. (2), (3), (7), and (9) in order to obtain: 

dY = (AC + ipAD) (xdt + dl t > T ) (10) 

where we have removed all the dependence in dS, just setting 8 — —dsC—ipdsC. 
In this case, we must ensure that either dY — 0, or dY has no definite sign. 
We will avoid the choice V> = —AC/AD, which leads back to the trivial case 
dY = 0. If dY > 0, and dY ^ 0, we will have to design a null portfolio whose 
value can only rise. In other words, this market shows arbitrage opportunities. 
Obviously, the reciprocal scenario, i.e. dY < and dY ^ 0, also lead to 
arbitrage, just building the Y = —Y portfolio. Clearly, the arbitrage possibilities 
other portfolios can be easily translated into Y terminology Therefore wc need 
to analyse the behaviour of dY, and thus we will find the constrains to \(t, S, a). 
Now it will be very convenient to recall the decomposition of dl t > T stated in 
Eq. (4), 

dY = (AC + ^AD) (( x + \l t <r) dt - XdG) , (11) 
and to inspect the properties of dG: 

[ t>T, 

dG= < dt- X- 1 t<T <t + dt, 
[ dt t >t + dt. 

It is clear that for t > r, dY reduces to dY = (AC + {pAD)\dt. But once the 
jump has happened, there is no financial reason to have a price that differs from 
the Black-Scholes price corresponding to a = o\,. And this is what we will get 
if we set \ — right after the change in the volatility. 

When the jump has not yet happened the differential dY reads: 



dY = (AC + tpAD) ( x dt + X(dt - dG)) , 
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with (dt — dG) > 0. Therefore we must choose x < for t < r. Collecting all 
this we obtain the following formula: 

X(t, S, a) = -fi(t, S) l a=aa = -Q(t, S) l t<T , 

where Q(t, S) is a strictly positive-definite bounded function depending on t and 
5. This function may depend, in a parametric way, on a a and at,, but as well 
on to and Sq. In fact, in general it may also depend on some parameters among 
those that characterise the contract specifications but, and this is a crucial point, 
never on all of them. We must have in mind that Eq. (9) must hold at least for 
another option D, different from C. Otherwise the market will not be complete. 

We have shown the mathematical properties that x(t, S, a) must fulfil al- 
though we have deepened very little in its financial interpretation. Let us intro- 
duce function ^(t, S, a), 

*(t,S,a) = (A-fi(t.S)) l t<r , 

in Eq. (11), 

dY = (AC + tpAD) {^dt - XdG) , 

and then evaluate the conditional expectation of dY, for a given value of S. Since 
E[dG] = 0, it is clear that E[dY\S] = E[V\S](AC + tpAD)dt. Thus *(t,S» 
measures the market price of the volatility risk, and it is exogenous to the market 
itself. It should be the financial agents who determine this function on the basis 
of their own appreciation of the actual risk. For instance, some authors [11] 
demand the absence of the so-call statistical arbitrage, that is _E[rfy|5 f ] = 0. 
This requirement implies that ^ — 0, i.e. that = A. 

5 Explicit solutions 

We can now solve Eq. (9) under appropriate conditions. For example, we shall 
begin assuming that O = A is constant, but not necessarily equal to A, 

d t C + ^<r 2 S 2 d 2 ss C -rC + rSd s C + \ACl t<T = 0. (12) 

This implies that the risk is felt uniform in time. We will also consider that the 
price of the option is constrained by the final condition: 

C(T,S,a;K) = $(S;K), 

which means that it will be a European-style option, where the price of the 
derivative in a fixed instant in the future, the maturity time, only depends on 
the actual value of the underlying at that moment and on some reference value, 
the strike, K. The function the payoff, will change for different kind of 
options within this same family. For instance, for the plain vanilla call we have: 

C(T, S, a; K) = m&x(S(T) - K, 0). 
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In addition, the mathematical nature of Eq. (9) demands that the solution 
satisfy two extra boundary conditions which, in this case, read 

C{t,0,a;K) = 0, and, lim C ^ S f' K ) = L 

Obviously, the same procedure can be used for other $ functions, such as the 
binary call where the payoff is 

C(T,S,a;K) = 1s(t)>k, 

where the boundary conditions to be fulfilled are 

C(t,0,a;K) = 0, and, lim C(t, S, a; K) = 1. 

S^oo 

Therefore, we will not specify a single function <!>, but we will treat all the 
suitable candidates at once. Moreover, we will use the term "Black-Scholes 
price , as a synonymous of the solution of the Black-Scholes equation for 

the given payoff, without further distinction. 

This will be the case when considering Eq. (12) for r < t, since then it 
reduces to the Black-Scholes model: 

dtC + \<J 2 h S 2 d 2 ss C ~rC + rSd s C = 0, 

whose solution is accordingly C(t,S,ab;K) = C BS (t, S, K). Nevertheless, 
we will show the main guidelines to solve it, because this will illustrate more 
sophisticated problems to come. The first step is to introduce two new variables, 
t* = T — t and x — log(5) + (r — of /2) (T — t), and to assume that C depends 
on its own arguments only through them: 

C(t, S, a h - K) - e-^-'V (t - t, log(S) + (r - ?pj (T - t);K 

This assumption implies the existence of a function of two variables V(t*,x; K) 
that obeys the following differential equation: 

d t *V =\<j 2 b dl x V. (13) 

Note that t* represents a reversion of the time arrow, that starts now at matu- 
rity. We have thus transformed our final condition into an initial one: 

V{0,x ;K) = ${e x °;K). 

This problem has a straightforward solution: 

/ + OO -, _(x-x ) 2 

dx $(e x °;K) e , (14) 

-oo ^2na;t* 
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and therefore, 



C(t, S, a b ; K) = e~ r ^ -*> / dx Q $(e Xo ; K) 



/+OC 
1 
-oo 



1 



(log(S) + (r- CT 2 /2 )(T-t)-x ) 2 



y/2na£(T-t) 

(15) 

which is the Black-Scholes price. When r > t the equation for C(t, S, a a \K) is 
a little more complex: 

d t C + ^a 2 S 2 d 2 ss C -rC + rSd s C + A (C BS (t, S, a b ; K) - C) = 0. 

The last term comes from the x^C contribution. The key point is to realize that 
in the expression for AC appears, not only C(t, S, <Jb;K), which we have found 
in Eq. (15), but also C(t, S 7 a a ; K), the unknown quantity. The procedure to 
follow is very similar to the one of the previous case. We will use again variable 
t*, and define £ as £ = log(S) + (r — cr 2 /2) (T — t). In fact x relates to £ through 
x = £ + (a 2 — a 2 ) t*/2, what will be useful in a forthcoming step. Now we 
assume again a particular dependence on the price of this new variables, 

C(t, S, a a ;K) = e -^ T ^Z (\og(S) + (r — \ (T — t),T — t] , 



where Z(t*,£;K) obeys the following equation, 

dt*Z = \<r 2 a d 2 a Z + Xe xe V ^ + ^d t *,t*;Kj , 
with the function V of Eq. (14); and the corresponding initial condition, 

After some algebra its solution reads 

r T 

C(t, S, a a ;K) = e - x( - T -^C BS {t, S, cr a ; K)+X J du e - x< - u -^C BS (t, S, a(u-t, T-t); K) 
where some short of "effective variance", d{t ai t b ), has been introduced: 

^(t a ,t 6 ) = < #° + °» ( * 6 ~* a) . (16) 

Note that a(0,T — t) = aj, and o(T — t,T — t) = a a . This behaviour can be 
used for compacting the solution. We can perform a typical integration by parts 
inside the integral sign and recover: 

C(t,S,a a ;K) = C BS (t,S,a b ;K) + J due - **"-*) d u C BS (t, S,a(u-t,T-t); K). 

(17) 
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The main benefit of the last expression is that it can be easily combined with 
Eq. (15), thus yielding: 



C(t, S, cr; K) = C BS (t, S 1 , a b ; K)+l t<T J du e'^^ d u C BS (t, S, a(u-t, T-t);K). 

(18) 

Note that this result is not sensitive to whether r is smaller than T, or not. 
Depending on the payofT function, the integral that appears in Eq. (18) can 
be computed, and analytic expressions for the option price can be obtained. 
In the most of the cases we have explored, however, the final formulas are 
cumbersome, thus providing little insight into the problem. We will present 
the simplest expression we have found by way of example. It corresponds to a 
vanilla call, in the special case that the discounted moneyness is equal to one, 
i.e. S = Ke~ r{ - T ~ l \ and that A, <j a and u b are such that they fulfil 

A=^i>0. 

8 

Then Eq. (18) reduces to 

C{t,Ke- T V-*\a\K) = C^^Ke-^^^^K) 



T-t 



2ir 

Up to this point we have reproduced the framework that corresponds to the 
problem stated by Herzel [7]. Our output agrees with his expression for t = 0, 
which is in fact the first option price given in Herzel's paper, although he does 
not number the equation. After that, he generalizes his formulation for any 
later instant of time, < s < T. Unfortunately, there is an erratum in his 
proposal. Thus, the limits in the definite integral in Eq. (4.28) should be and 
T — s, instead of s and T. Or, in an equivalent way, t should be replaced by 
t — s, keeping the rest unchanged, including dt. 

More general solutions can be obtained using the same approach, with little 
extra effort. We can consider, for instance, the case of a % depending on all the 
involved time magnitudes: 

x = -v(t;to,T)i t<T , 

with n(t; to,T) > 0. The solution for t = to is simply 6 : 
C(t ,So,v,K) = C BS (t ,So,<T b ;K) 

+ l t0<T due- J "o dtTl{t]to ' T) d u C BS {to,So,a{u-to,T-t );K). 

Jt 

(19) 



'Given that, up to this moment, the notation did not induce to misunderstanding, we had 
not stressed the difference between to, the actual time in which the options is evaluated, and 
t, a generic instant of time, tg < t <T. 
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We have thus obtain a broad set of valid prices with no clear financial inter- 
pretation. A plausible requirement that may help us to discard candidates is 
to demand that the final solution only depends on T — to- This forces that 
r](t;to,T) = f(t — t ;T — to)- In Appendix A we find that, if we follow an 
heuristic approach that tries (in vain) to cancel out all the risk, the option price 
is 

C(to, So, a; K) = C BS (to,So,cr b ;K) 

r T 14- e — A(it— 1 ) 
+ l t<T du d u C BS (to,S ,a(u-t ,T-to);K). 

Jto 1 

(20) 

We will disregard at this moment the possible interpretation of this solution, 
and concentrate in its validity instead. It is straightforward to check that if we 
replace 

e -A(t-to) 

„(t; to, T) = V (t t ) = A i + e _ A( ^ Q) , (21) 
in Eq. (19), we will recover Eq. (20). 

6 Numerical computation 

Now it is time to analyse and compare the different solutions we have found, 
and eventually to represent some of them. This task is easier if we express our 
results in the form of the expected value of the discounted payoff, under some 
appropriate probability density function: 

C(t Q ,So,* ;K) =EQ[e- r ( T - t °^(S(T);K)]. (22) 

Thus we will be able to translate the functional form of \ m to an equivalent 
model. Since all the possible candidates to be the fair price collapse to the 
Black-Scholes solution if the jump takes place, we will assume that t < t from 
now on. 

Let us begin with 77 = A. In that case the final price S(T) in Eq. (22) can 
be expressed in the following terms 

S(T) = S oe (^/ 2 )(--M+-aW(r-t ) lt< _ 

+ S(f)e( r - a2b/2 ) {T - f)+ab [ W(T - to) - W{f - to) h t>f , 

where, as usual, fi has been replaced by r, and W(t — t ), a new Brownian 
motion with zero mean and variance equal to t — to, has been introduced. The 
jump process f follows also a different exponential law: 



P(to<r <t) = l-e- X{t - tQ \ 
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Figure 1: Option pricing in terms of the present moneyness, Sq/K for two 
different maturities: (a) T — to = 0.25 years, and (b) T — to = 0.05 years. The 
depicted results correspond to a vanilla call, &(S;K) = max(5' — K,0). The 
numerical value of the involved parameters are r = 5%, a a = 10%, Ob = 20%, 
A -1 = 0.1 years, and K = 100, in suitable currency units. The "constant" 
line correspond to setting rj = A, whereas by "heuristic" we mean the choice in 
Eq. (21), which was somewhat inspired by an heuristic approach. Note that the 
first price is more similar to the plain Black-Scholcs price with a = at, arid that 
conversely the second method leads to a price closer to the Black-Scholes one 
for a — <r a . The discrepancy is reduced as the maturity time approaches. All 
the plots were obtained using Monte Carlo techniques over 100 000 replicas. 
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Thus when A ^ A we will consider in practice that the model is not accurate in 
the forecast of the actual mean transition time, and that it should be replaced 
by another value. It is not our intention to add superfluous complexity to the 
study of the several proposals we have done. Therefore we will concentrate in 
the case that the original A (and r) is used. 

Summing up, we can actually compute the price of the option under the 
previous assumptions with the following procedure. We choose a value for r that 
follows the proper probability density. If that time is bigger than the maturity 
we need to generate a zero mean Gaussian variable, with T — to variance. Thus 
S(T) will be governed only by a a . Conversely, if r < T we will need two 
independent zero mean Gaussian variables, with variances equals to t — to and 
T — t. In this case a a governs the behaviour of the equivalent stock price up 
to S(t), and a\> does thereafter. The averaged value of the discounted payoffs 
will lead to the correct estimation of the desired magnitude. This set-up is 
the appropriate to perform Monte Carlo numerical simulation, as it is shown in 
Fig. 1. 

The choice in Eq. (21) for r\ is equivalent in law to consider that 
S(T) = Soe( r -^ /2 )^ ta)+ ^ W{T - to h t<T 

+ l5( T ) e ('--^/2)(T-r)+<7 [w : (T-to)-W : (r-to)] lt>T 
+ l5( T ) e ('--^/2)(T-r)+<7 6 [w : (T-to)-W : (T-to)] lt>r) 

with the same specifications for r and W(t — to) as in the past scenario. Note 
that this particular choice for the market assessment of risk leads to the follow- 
ing alternative strategy. We choose again a value for r, using an exponential 
probability density function of mean 1/A. If that time is bigger than the matu- 
rity, nothing changes with respect to the previous example, and a a drives the 
evolution of the underlying all the time. But when t < T the picture changes. 
The price S(T) is the arithmetic mean of the two possible paths: one in which 
the volatility is a\, right after the jump, and the other that considers that a a 
remains unchanged. This is the reason to dissociate the premise that assumes 
the existence of a distinctive time value r, and the innovation that it carries. We 
can be sure that forthcoming news may affect the market and, at the same time, 
only guess about the final effect. Thus Eq. (21) leads to a more conservative 
risk analysis, in the sense that this price is near the Black-Scholes value corre- 
sponding to a a , whereas the 77 = A choice anticipates more intensely the future 
change in the volatility. This explains the behaviour of the different call prices 
observed in Fig. 1. Also in this figure, but specially in Fig. 2, we can check that 
the two prices converge to the no-jump solution as the maturity horizon comes 
closer. 

Finally, we want to point out that we can obtain solutions in all the range 
of confidence levels for the model. Thus, if we set q equals to the probability 
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0.00 0.05 0.10 0.15 0.20 0.25 
Time to maturity 



Figure 2: Option pricing in terms of time to maturity when the present mon- 
eyness is one. We consider an at-the-money call in the same set-up described 
in Fig. 1. The plot shows how both prices change when the time to expiration 
shortens. Clearly the price that fully ignores the risk of a change in the volatility 
becomes more accurate as the probability of a jump reduces. 
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that the volatility actually varies when the jump takes place, the function 77 that 
express this risk evaluation is 

e -A(t-to) 

V (t;t ,T) = V (t-t ) = X {i _ q)/q + e _ x{t _ toy 

The only constrain is that we cannot neglect all the risk just setting q = 0. 
Therefore < q < 1 leads to a valid price, although in absence of further 
information q — 1/2 and q = 1 seems to be the only privileged values. 

7 Conclusions 

We have revisited the framework stated by Herzel, in which the dynamics of one 
asset S is driven by a log-normal diffusion equation with an stochastic volatility 
parameter a. The volatility of this stock may jump at a random time r from a 
fixed initial value a a to another fixed final value ct&. And no more than one of 
such jumps is allowed. This event can model, for instance, the future publication 
of crucial information related to this specific market. 

We have introduced a procedure for obtaining fair option prices, different 
from the used one in Herzel's manuscript. There, the author exploits inten- 
sively probability arguments for finding the necessary and sufficient conditions 
that the model must fulfil to be complete and arbitrage-free. Thus he derive 
the equivalent martingale measure Q. We have employed another technique 
of broad use in this field. We have determine the partial derivative equations 
that, according to the Ito convention, the option price must fulfil. We have 
shown that the use of a secondary option completes the market. After that, 
we have demanded that the market has no arbitrage and we have found the 
exogenous function that measures the market price of the volatility risk. We 
have explored the output for several choices of this function and, incidentally, 
we have amended some of the results presented in the original reference, where 
the risk premium was null. 

In fact, one of the biggest benefits of our approach is when considering more 
sophisticated prescriptions for the market price of volatility risk. We have not 
only obtained closed formulas in such a cases, but we have also been able to 
interpret the financial meaning of them. We have seen how a choice for the 
volatility risk price can be translated into a lack of confidence in the model 
premises. For instance, a constant risk price, other than zero, plays the same 
role of a redefinition in the mean transition time of the jump process. 

In particular, we have studied with some detail a solution that can be under- 
stood as the response of a suspicious market maker, who admits the possibility 
that volatility stays in the same level, although the jump (that is, the announce- 
ment) has taken place. We have also presented plots with actual examples of 
these solutions, computed using Monte Carlo numerical techniques. 
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A Appendix 

In this appendix we present an heuristic approach trying to remove the risk 
associated to the volatility change in a portfolio without secondary calls. It 
it obvious that we will fail in this task by construction, because this is only 
feasible if the jump time is deterministic. But in this case, Eq. (2) and the 
expressions derived from it are not longer correct. Therefore, we will assume 
that t is a stochastic magnitude, and the solution we obtain will be conditioned 
to its value. Let see how it works just recovering Eq. (10), and setting 'tp = 0: 

dY = AC ( X dt + dl t > T ). 

We are now demanding that dY = 0. This constrain leads to: 

X=~f t U>r = -8(t-r), 

and to the corresponding equation is then: 

d t C + ^a 2 S 2 d 2 ss C -rC + rSd s C + ACS(t - r) = 0. (A.l) 

Note that this equation corresponds to Eq. (6) with S = dsC and ip — 0, as we 
have just stated. We will search for a solution of the form 

C(t, S, a; K) = e^^U (t - t, log(S) + [r - y ) (T - t); K 

based upon the two variable function U(t*,x; K), which must fulfil the following 
partial equation: 

d t ,U-±* 2 d 2 xx U = AU6(T-T-t*). 
We therefore consider the Fourier-Laplace transform of x; K), 

r-\-oo r-\-oo 

U(s,lu;K)= dt*e- st * dxe iuJX U(t*,x;K), 

JO J-OO 
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that follows the simpler equation: 

sU-U + \f 2 u?U - AU T - T e~ s{T ~ T) lr<T = 0, 

where the tilde stands for the Fourier transform of the corresponding object. 
Thus Uo(iv; K) = U(t* = 0, u, a; K), and it does not depend on a. On the other 
hand, AUt- t = U(T — t,u, <Jb', K) — U(T — r, u, <r a ;K). Now we can isolate all 
the explicit dependence on the Laplace variable s, 

U(s, u>, a; K) = \ f U + AU T - T e-^ T -h T<T } , 

and perform an inverse transformation, 

U(t*,u;,v;K) = U e- M */ 2 + AU T _ T e-^ 2 ( e - T+ ^/ 2 l T _ t , <T < T . (A.2) 

Notoriously, the second term only gives contribution when the jump is comprised 
between t, and the maturity, T. When t* < T — r, i.e. r < t and a = a^, 
Eq. (A.2) reduces to, 

U(t*,uj,a b ;K) = U^> H ' I 2 - U BS (t*, u, <J b ; K), 

which leads to 

C(t,S,a b ;K) = C BS (t,S,a b ;K), 

a riskless price. 

When t* > T — r, that is to say, when t < r and a — <r a , but r > T, the 
main equation also takes a simple form, 

U(t*,cj,a a ;K) = U*e-°^ H ' ' 2 = U BS (t* ,tj,a a ; K). 

In this case, since the change in volatility occurs after the maturity of the 
contract, the price reduces to a plain Black-Scholes model without any jump of 
volatility, 

C(t, S, <i a ; K\T < r) = C BS (t, S, a a ; K), 

Thus, this scenario has again no risk associate with it. 

Finally, when t* > T — r and r < T, all the terms contribute to a more 
complex expression, 

U(t*,uj,a a ;K) = U^r^ 1 ' I 2 + AUr-re-^ 2 ^' - T+T ^ 2 . (A.3) 

Recall that AUt-t = U(T — T,w,a b ;K) — U(T — T 7 uj 7 a a ;K) is a term that 
counts only for the variation in U due to the change in the volatility, when it 
takes place. Thus U(T - r, uj, a h ; K) = U BS (T - t, uj, a h ; K). The other term 
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can be obtained by self-consistency. We will start from Eq. (A. 3) and take a 
limit: 



U(T - T,Lu,<r a ;K) 



which leads to 



lim U(t*,u},a a ;K) 

t* — >T — T 

lim U e-^ 2t '/ 2 + AU T -re- a ^ 2it '- T+T)/2 , 



U(T-T,LU,a a ;K) = - [U BS (T - r, u, a a ; K) + U BS (T - r, u>, a b ; K) 
Now we will introduce this result back into Eq. (A. 3), and obtain so 



U(t*,uj,a a ,K) = -U 



<J a u 2 t' /2 _|_ e -5 2 {t* -T+T,t*)u 2 t' /2 



where a(t a , t b ) is the same function which we have previously defined in Eq. (16). 
Therefore 

C(t, S,a a ;K\t<T<T) = ± [C BS (t, S, a a ; K) + C BS (t, S, a(r - t,T - t);K)} , 

Finally, in order to obtain a expression for t < t that does not depend on 
future information, we will compute the expected value of the previous condi- 
tioned solutions: 

C{t,S,a a ;K) = E[C{t, S,a a ; K\r = u)} 

= \J du[C BS (t,S 7 <j a -K) + C BS (t 7 S 7 <j(u-t,T-t);K) 



-X(u-t) 



A 



+ 00 



duC BS {t,S, a a ;K),e 



-\(u-t) 



an expression that reduces to 
C(t,S,a a ;K) = C BS (t, S,a b ; K) 
Then the complete result is 



T i , -X(u-t) 

BSuc^..Tr\, / du—— d u C BS (t,S,a(u-t,T-t);K). 



C(t, S, a; K) = C BS {t, S, a b ; K)+l t<T J du L±£_^! - d u C BS {t, S, a(u-t, T-t);K), 
which does not fulfil Eq. (A.l), but it is still a valid solution. 
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